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FOREWORD 


This  report  describes  the  results  of  some  preliminary  theoretical  work  in  the  area  of 
quantum  mechanical  computing  which  has  been  performed  by  the  Strategic  Systems  (K) 
Department  at  the  Naval  Surface  Warfare  Center  (NAVSWC)  as  part  of  a  molecular 
computing  focused  technology  initiative. 

This  report  has  been  reviewed  and  approved  by  Ted  Sims,  Space  Sciences  Branch 
Head,  and  James  L.  Sloop,  Space  and  Surface  Systems  Division  Head. 


Approved  by 

R.  L.  SCHMIDT,  Head 
Strategic  Systems  Department 
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INTRODUCTION 


Despite  the  obvious  practical  considerations  (e.g.,  stability,  controlability),  certain 
quantum  mechanical  systems  seem  to  naturally  lend  themselves  in  a  theoretical  sense  to  the 
task  of  performing  computations.  The  purpose  of  this  report  is  to  describe  one  such 
idealized  system-the  quantum  harmonic  computer.  As  its  name  might  suggest,  this 
theoretical  device  employs  the  well  known  energy  characteristics  of  quantum  mechanical 
oscillators;  the  associated  creation  and  annihilation  operators;  and  the  quantum  mechanical 
axioms  of  state  preparation  and  observability  to  perform  computations.  It  is  demonstrated 
that  programs  can  be  written  for  this  device  in  terms  of  quantum  mechanical  observables 
and  creation  and  annihilation  operators  which  will  algorithmically  manipulate  oscillator 
energy  states  to  perform  the  desired  calculations,  the  results  of  which  are  eigenvalues  of  a 
well  defined  system  observable.  By  definition,  these  programs  are  equivalent  to  Turing 
machines,  so  that  anything  that  is  Turing  computable  is  also  computable  with  this  device. 

It  is  well  known  that  many  physical  systems  (e.g.,  photon  gas,  crystals,  phonons,  diatomic 
molecules)  can  be  approximated  by  collections  of  quantum  harmonic  oscillators.  Although 
the  development  herein  may  suggest  such  implementations,  we  emphasize  that  our  intent 
has  been  theoretical  in  nature  rather  than  practical. 


QUANTUM  MECHANICAL  HARMONIC  OSCILLATOR 


In  order  to  provide  a  more  complete  foundation  for  the  discussions  which  follow,  we 
provide  in  this  section  a  brief  sketch  of  the  one-and  m-dimensional  quantum  mechanical 
harmonic  oscillators.  Although  we  adopt  the  Dirac  notation,  our  development  is  standard 
and  can  be  referenced  in  more  detail  in  most  basic  quantum  mechanics  texts  [e.g., 
References  1,  2]. 

The  ideal  classical  harmonic  oscillator  is  a  particle  of  mass  m  constrained  to  move 
along  an  axis  and  subject  to  a  restoring  force  proportional  to  the  displacement  from  the 
force  center.  The  corresponding  quantum  mechanical  oscillator  is  a  particle  of  mass  m  in 
one  dimension  with  energies  generated  by  the  Hamiltonian  operator 

H=— (p2+m2(j32q2),  (2-1) 

2m 
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where  w  is  the  natural  angular  frequency  of  oscillation,  and  q  and  p  are  the  position  and 
momentum  operators,  respectively,  related  by  the  commutator 

[&p]=ik  (2-2) 


Here  h  is  Planck’s  constant  divided  by  2n. 

The  energy  eigenvalue  problem  for  this  system  is  given  by  the  time-independent 
Schrodinger  equation 

H|E>=E|E>,  (2-3) 

where  E  is  the  energy  eigenvalue  associated  with  the  system  state  described  by  the  "ket" 
vector  |E>. 


Using  this  we  may  now  construct  the  vector  space  of  the  dynamical  states  of  the 


system.  Let 

d=(2hm(ji)~l,2(p-im(dq ) 

(2.4) 

and 

a+=(2>iW(o)"1/2(p+im(og). 

(2.5) 

It  is  readily  seen  that  the  operator  a  and  a*  are  Hermitean  conjugates  and  that  application 
of  (2.2)  provides  the  commutator  relation 

[641= 1.  (2-6) 

Furthermore,  it  is  easily  shown  that  the  Hamiltonian  operator  of  (2.1)  can  be  written  in 
terms  of  6  and  a*  as 

H=(N+l/2)h<o,  (2-7) 


where 

N=a+fl. 


(2.8) 


Thus,  the  eigenvalue  problem  of  (2.3)  is  equivalent  to  constructing  the  eigenvectors  of  the 
operator  N.  It  should  be  noted  that  N  is  Hermitean  and  a  system  observable. 

For  the  sake  of  brevity,  we  omit  its  detailed  development  and  simply  state  the  well 
known  result  that  the  spectrum  of  eigenvalues  of  N  is  the  set  of  non-negative  integers  {0, 

1,  2,  ...,  n,  ...}  from  which  we  may  form  the  orthonormal  state  vectors 

|0>,  1 1>,  |2>,...,  |  n>,...,  where 

N|n>=n|/i>.  (2-9) 

For  the  same  reason  as  above,  we  also  state  without  proof  the  following  well-known 
properties  of  the  operators  6*  and  6: 


a*\n>={n+\)ll2\n+\>  (2.10) 
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(2.11) 

(2.12) 


a|/i>=(«)1/2|n-l> 
a|0>=0 
and 

|n>=(n!)‘1/2(a+)',|0>,  (2-13) 

where  |0>  is  called  the  ground  state  of  the  system. 

Upon  application  of  (2.7)  and  (2.9)  to  (2.3),  it  is  easily  verified  that 

H|7t>=(n+l/2)ho)  \n>,  (2-14) 

i.e.,  the  energy  states  of  the  oscillator  are  quantized  into  discreet  multiples  of  hw  with 
ground  state  energy  l/2hw. 


Consider  now  the  m-dimensional  system  of  m  distinguishable,  noninteracting 
oscillators  with  system  Hamiltonian 


(2.15) 


where  each  ft,  in  the  summation  has  the  form  of  (2.7)  with  all  associated  operators  having 
the  appropriate  i  subscript.  If  |  nj>  denotes  the  eigenvectors  of  H;,  then  the  tensor  product 
of  eigenvectors 


(2.16) 


forms  a  complete  orthonormal  set  Sm  of  eigenvectors  of  H  in  (2.15)  so  that 

m 

H|»i«2-nm>=£  (fti+mhu)  |«1n2...nm> 

i=l 

and 

m  m 

N| »i.n2  •»m>=E NiK*r..nm>=E n.\nxn2...nm>. 

M  i=l 


(2.17) 


(2.18) 


Analogous  to  (2.6),  it  can  be  shown  that  the  operators  at  and  d-  obey  the  following 
commutation  rules: 
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[<*>;]  =0 

[«i\«/]=0  ►  (ij  =  0,1, .•,»»)  (2.19) 


It  is  convenient  at  this  point  to  introduce  several  new  operators  which  will  prove 
useful  in  the  following  sections.  Define  the  new  operator  (3 .  as 


P,-h>=| 

[  (nj)"u2|ni>,  n,*0 
{  0  ,  0 

(2.20) 

and 

a  i s  P  fli 

(2.21) 

A 

A  A  A 

araiPj 

(2.22) 

so  that 

a  + 

A 

+ 

gT 

ii 

A 

g~ 

(2.23) 

and 

&i\ni> 

||nrl>,  nt*  0 

1  0  ,  «.=  0 

(2.24) 

It  is  easily  verified  that  &  .  and  a)  obey  the  following  commutation  relation  for 
1  <.i,  j<m: 


1,  i=j,  nt  =  0 
0,  otherwise 


(2.25) 


We  also  define  the  identity  operator  fj  as 


so  that  for  any  operator  y  t  , 

Y.WrY,  • 


(2.26) 

(2.27) 
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QUANTUM  HARMONIC  COMPUTER:  DEFINITIONS 


In  this  section  we  develop  a  set  of  definitions  that  will  provide  for  a  precise  description 
of  the  quantum  harmonic  computer.  Let  Sm  be  as  defined  in  the  previous  section  and 
QmcSm>  where 

Qm={\nyn2...nm>  \nie{0,l},3zizm;3zn2zm).  (3.1) 


Each  eigenvector  in  Qm  is  called  an  instantaneous  description  (ID).  The  quantum  numbers 
nx  and  n2  in  each  ID  are  called  the  state  and  pointer,  respectively.  Also  define  the  following 
sets  of  operators: 


where 


and 


m={ap& ]\UUm), 

(3.2) 

={N,,  Nj,  ...Na,  T) 

(3.3) 

m 

t=EN,. 

(3.4) 

i=3 

R  =£fmU  T|  . 

(3.5) 

We  note  that  f  is  Hermitean  and  is  a  system  observable. 

An  Rm-quadruple  is  a  4-tuple  of  one  of  the  following  three  types: 

Ak  l 
n.n„A„  A 
1  »2  n2  r 

A  +  .  Z 

ni\a2AV 


l 


(3.6) 

(3.7) 

(3.8) 


where  is  the  quantum  number  describing  the  eigenstate  of  the  n2  oscillator;  k  and  /  are 
non-negative  integers;  A;eSfm;  and 
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{Ai  Ai  ...  i4j.  ,  k* 0 

'ktimes  (3-9) 

/,  ,  k=0. 

Each  i?m-quadruple  may  be  interpreted  as  being  the  instruction:  "if  the  state  is  nx  and  the 
eigenstate  of  the  n2th  oscillator  is  for  some  ID,  then  apply  the  following  operators  to  the 
ID." 


A  program  Pm  is  a  finite  nonempty  set  of  Pm-quadruples,  no  two  of  which  have  the 
same  state  quantum  number.  Let  p-n^n^  represent  any  Pm-quadruple,  where  x  is  any  of 

the  operator  pairs  in  (3.6)  -  (3.8),  and  X,  Y  e  Qm.  A  program  Pm  induces  a  basic  Pm 

transition  X^-+Y  from  X  to  Y  if  there  is  a  pePm  such  that  N.X  =  n.X;  N  X  =  n  X:  and 

Y  =  f  X.  A  finite  sequence  Xj— 7  of  such  transitions  is  called  a  /^-computation 

if  there  is  no  p"ePm  or  WeQm  with  Z^—+W.  In  this  case  the  eigenvector  Z  is  called  the 
resultant.  The  eigenvector  X  is  called  the  initiator.  It  will  always  be  assumed  here  that  for 
every  initiator,  the  pointer  quantum  number  n2  =  3  and  the  state  quantum  number  ^  =  1. 

We  may  now  define  the  m-dimensional  quantum  harmonic  computer  rm  as  the 

4-tuple 

«  (MO) 


Informally,  rm  can  be  thought  of  as  an  m-dimensional  quantum  mechanical  harmonic 
oscillator  which  serves  as  a  memory/storage  medium  for  any  ID  in  Qm.  This  oscillator 
interfaces  with  a  "device"  which  imposes  upon  it  energy  eigenstate  transitions  via  the 
application  of  a  sequence  of  operators  from  Rm  as  algorithmically  dictated  by  Pm. 


NUMERICAL  COMPUTATION  USING  rm:  SIMPLE  EXAMPLES 


It  is  obvious  that  there  is  a  one-to-one  correspondence  between  Pm  programs  and  Turing 
machines  [3,  4].  Hence,  as  long  as  we  view  m  as  being  potentially  infinite  (i.e.,  adding  extra 
oscillators  as  needed),  anything  Turing  computable  is  also  Pm  computable.  In  this  section 
we  provide  examples  of  Pm  program  which  compute  several  non-negative  integer  arithmetic 
functions  and  illustrate  the  associated  Pm  computations  using  simple  initiators. 

In  order  to  perform  these  computations,  we  make  use  of  the  following  symbolic 
representation  for  a  non-negative  integer  J: 


6 


NAVSWC  TR  91-205 


/=  11...1, 

~  J+ 1  times 


(4.1) 


(thus  5'  ■  111111).  Also,  we  may  assume  without  violating  any  quantum  mechanical 
principles  that  initiators  may  be  prepared  with  the  required  state  and  pointer  quantum 
numbers,  as  well  as  any  combination  of  quantum  numbers  n^ {0,1}  for  3<i<m.  Initiators 
may  therefore  be  prepared  which  can  contain  within  the  energy  states  of  oscillators  3 
through  m  symbolic  representations  for  non-negative  integers  of  the  form  (4.1).  We  adopt 
the  convention  that  if  /  and  K  are  two  non-negative  integers  which  are  to  be  symbolically 
represented  within  an  initiator  with  J  appearing  first,  then  the  initiator  will  be  prepared  with 
nj  =  1  for  3  <.  i  <.  J+3,  J+5  <  i  <.  J+K+5,  and  n4  =  0  elsewhere  for  i  >  2  (thus  m  >. 
J+K+5).  An  initiator  prepared  in  this  manner  is  said  to  be  properly  prepared. 

Each  of  the  sample  Pm  programs  described  below  are  designed  to  produce  a  resultant 
from  a  properly  prepared  initiator  via  the  sequential  application  of  operators  in  Rm.  The 
computed  solution  is  the  eigenvalue  of  T  when  applied  to  the  resultant. 

Example  1.  The  set  A  is  a  Pm  program  which  computes  J+K  for  J,K  non-negative  integers, 
where 


A^lld^d^0,  lOdjdi,  21d2(a1)°,  Sld^d^0}. 


For  the  sake  of  clarity  we  have  used  (dj)0  in  A,  but  note  that  by  (3.9)  and  (2.27)  it  need  not 
be  made  explicit.  Effectively,  this  Pm  program  produces  a  resultant  in  which  oscillators  3  and 
/+5  have  transitioned  to  their  ground  states  and  (neglecting  the  state  and  pointer)  J+K 
oscillators  remain  in  their  first  excited  state.  Thus  the  computed  solution  J+K  is  the 
eigenvalue  of  f  when  applied  to  the  resultant. 


In  order  to  illustrate  this  let/=2  and  K- 
Ps  computation  is: 

d3(d1y 


-1  and  choose  m=J+K+5=8.  The  associated 


|13111011>- 


►|13011011: 


A  +  A  + 
«2«1 


12401101 1> 


|25011011> 


a2+(«i)° 


A  +  A  -• 

a2al 


•|26011011>— — -*|37011011> 


d7(  d/ 


1 3701 1001 >, 

so  that 
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T  |3701 1001>=3 1 3701 1001>. 


As  an  aside,  we  note  that  in  general  the  composite  operator aJ+Ja3(&2)y+2(&i)2will  always 
produce  the  desired  resultant  for  J+K  for  any  properly  prepared  initiator. 


Example  2.  The  set  B  is  a  Pm  program  which  computes  J-K  for  J>K  where  J  and  K  are 
non-negative  integers: 


8=(n%(«|)‘.  20 314j(4,)%  3042«;,  414.(4,)",  404,4;,  51424;, 


«  W.  6042«;,  714,4;,  704, *(4,*)2,  814,(4/,  804*(4,)7,  904*(4,)“,  914,(4,)°}  . 


Although  not  necessary,  we  have  again  made  (c^)0  explicit.  Here  B  produces  a  resultant 
in  which  oscillators  3  through  K+3  and  7+5  through  J+K+5  have  transitioned  to  their 
ground  states  and  J-K  oscillators  remain  in  their  first  excited  state.  Thus  the  eigenvalue  of 
T  when  applied  to  the  resultant  is  J-K.  B  also  needs  a  ground  state  oscillator  after  the 
representation  of  K  in  the  initiator.  Thus  we  choose  m=J+K+6. 


As  an  illustration,  let  7=2  and  K=  1,  as  before,  so  that  m=9.  P9  produces  the 
following  computation: 


| 1311101 10> 


|130110110: 


<x2a, 


|2401 101 10> 


a2+  («i)e 


|250110110: 


&KK )° 


A  +  A  + 

a2al 


1 2601 101 10> - ►  1 3701 101 10> 


a2(«i)° 


|380110110> 


|390110110: 


a2al 


|4801 101 10> 


«8(«!)0 


|4801 10100>— ' |5701 10100> 


a2ai 

-J— ►  |6601 10100>^-4- 


|750110100>i^|840110100>-^lL|830110100>-^£. 
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A  /  A  \  O  A  A  + 

|47CX)10100>  —  -  |470010000>-^^*- 1560010000>, 

so  that 

T|560010000>=1 1560010000>. 

Again  we  note  that  the  composite  operator 

^/+iC+5•••^y+5®JC+3•••®3(“2)J,  (“l) 

will  always  produce  the  desired  resultant  for  J-K,  J>K,  for  any  properly  prepared  initiator. 

r4  LOGIC  GATES 

Because  of  the  obvious  equivalence  between  Turing  machines  and  Pm  programs,  we  can 
construct  Pm-programs  which  produce  logic  function  computations.  In  this  section  we  prove 
a  series  of  theorems  which  state  that  specific  4-dimensional  quantum  harmonic  computers 
can  serve  as  basic  logic  gates.  Let  us  assume  that  the  energy  states  of  the  third  and  fourth 
oscillators  represent  truth  values  for  propositions  r  and  s,  respectively,  where  n—O  means 
"false"  and  n,=  l  means  "true"  for  ie  {3,  4}.  The  truth  value  obtained  from  a  r4  logic  gate 
computation  is  the  eigenvalue  of  T  when  applied  to  the  associated  resultant. 

Theorem  1.  (Q4,  R4,  V,  X)  is  an  OR  gate,  where 

v={iia2x,  iod2+(«;)2, 20(d2rd;,  21*^;}. 

Proof: 
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£  +  A  +  A  A  + 

cutti  a, a, 

1 1311> - ►  |2411>-±-U  |3410>;f  |3410>=1  |3410>. 


&2&1  (a/al 

1310> - ►|2410> — - +  |3410>;f|3410>=l|3410>. 


i  a2(ai)2 

1301>— — — — >  |3401> 


;f|3401>=l|3401: 


A  +  /  A  +\2 

,  a2(cc.) 

1 1300>— — — ►  |3400> 


;T|3400>=0|3400>. 


Theorem  2.  (Q4,  R4,  A  ,  A)  is  an  AND  gate,  where 


Proof: 


A=jlla2aJ,10a2(ai)2,21an2(cti)2,20a2(di)3,31d>Ijai,30(aJl2)oaJ,51dB2dj|. 


cL(ttt)2 

1311> - *  |2411>  ►  |4410>;f  |4410>=1  |4410>. 


&2&1  a0(a|)3  a  a 

1 1310> - ►  |2410> — - — - — ►  |5310> — ^ ^  |4300>;f  |4300>=0|4300>. 


,  aU&l)2  &X 

1 1301> - ►  |3401>^-U  |4400>;f|4400>=0|4400>. 


1 1300> - ►  |3400> — — — -*  |4400>;f  |4400>=0 14400>. 
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Theorem  3.  (Q4,  R4,  V  ,X)  is  an  XOR  gate,  where 


^|llA2«i»10«2(^i)2>21®  j^)0  (&i)3,31(d  Bj)°  (&i)2,30(a  ^)0  (&i)2,40d2(di)2,61a  ^djj. 


Proof: 


1 131  |2411>  g4(a--')->  |4410>  g2(g--->  |6310>-^^  |5300>;f  |5300>=0 15300>. 


. <«4)°(*I)3 


1310>— - — ►  |2410> — — - |5410>;T|5410>=1 15410>. 


1 1301>a2(CCl)  >  |3401>(tt4)  |5401>;f  |5401>=1 15401>. 


,  d2(a;)2  (d^Cdi)2 

1 1300>— — — ►  |3400>  V  »  |5400>;T|5400>=0 15400>. 


Theorem  4.  (Q4,  R4,  ~V,  X)  is  a  NOR  gate,  where 


~v=(iid^d;,  iod2+(di)3, 20d2+d;,  sid^fdD2,  loi&^n&D2, 4 la^a;,  40d;d;}. 


Proof: 


A  A  +  ^  +  A  + 

a3al  ,  «2al 

*3  1 .  -i  LI 


d4(d;)2 


131 1> — = — ►  |2301>— - — *-|3401>— — — ►|5400>;7’|5400>=0|5400>. 


,  Mi  Mi  (dj°(di)2 

1 1310>— — ►  |2300>— — j 3400>  — — M  |5400>;T|5400>=0 15400>. 


a2(ai)3  a4ai 

1301>— - - ►  |4401>  >  |5400> 


;f|5400>=0|5400>. 
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1 1300>  g2(<Xl)3,  |4400>-^  |5401> 


;f|5401>=l|5401>. 


Theorem  5.  (Q4,  R4,  ~A,  X)  is  a  NAND  gate,  where 


~A={ll«2+a;,  10a2(&^21a^a^0(a^)o(«^30a2a;,41d^(a^51(a^)oa:,50a;a;}. 


Proof: 


1 131 1> 


a2«!  _  Mi  .  M 


4ai  i  cc,a,  <L(a,r 

241 1> — *-►  |3410> - -*  |4310>— - ►  |6300>;f  |6300>=0 16300>. 


,  &2&,  (aJ°(M4 

1 1310> - ►  |2410>— — - U  |6410>;f  |6410>=1 16410>. 


1 1301>— 2(Ctl)  »  |5401>(tt4)  |6401>;f  |6401>=1 16401>. 


1 1300> - ►  |5400> - +  |6401>;f  |6401>= 1 16401>. 


Theorem  6.  ( Q4 ,  R4,  ~V  ,  X)  is  a  NXOR  gate,  where 


~y=  {ll&^I.lOd^d^Od^;,  30(  aji&l)2,  31  41d|l2a1+,  40a;&:} . 


Proof: 


1 1311> — |2301>-^U  |3401>(tt4)  (0Cl)  ►  |5401>;f|5401>=l  |5401>. 


12 


NAVSWC  TR  91-205 


1310> 


Mi  ,  Mi  (CLJ°(&T)2 

_i_U  |2300>— — — >  |3400>— - 1 5400>;7’| 5400> =0 15400> 


A  +  /  a  +\3 

1 1301  >  —  —  »  |4401>— 1*. |5400>;f|5400>=0|5400>. 


a4al 


oti(tti)3  aiai 

1 1300>— — — ►  |4400>— — — *■  |5401>;7’|5401>=1 15401>. 
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